Evolution of magnetic and crystal structures in the multiferroic FeTe 2 05Br 
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Neutron diffraction and nuclear quadrupole resonance (NQR) measurements were employed to 
investigate magnetic order in the non-ferroelectric phase preceding the low-temperature multiferroic 
state in FeT^OsBr. Refinement of the neutron diffraction data and simulations of 79 ' 81 Br NQR 
spectra reveal that the incommensurate magnetic ordering in the non-ferroelectric state comprises 
amplitude-modulated magnetic moments, similarly as in the multiferroic state. The two ordered 
states differ in the orientation of the magnetic moments and phase shifts between modulation waves. 
Surprisingly, all symmetry restrictions for the electric polarization are absent in both states. The 
different ferroelectric responses of the two states are thus argued to arise from the differences in the 
phase shifts between certain modulation waves, which cancel out in the non-ferrolectric state. 

PACS numbers: 75.85. +t; 76.60.-k; 75.25.-j; 75.30.Kz; 71.15.Mb 



I. INTRODUCTION 

The discovery of magnetically induced electric 
polarizatioif^ revealed a new aspect and a great applica- 
tion potential of geometrically frustrated spin systemsP^ 
Frustration often leads to complex incommensurate (IC) 
magnetic structures, which can break the inversion sym- 
metry and thus overcome a fundamental restriction for 
the macroscopic electric polarization. Depending on the 
magnetic order, the magnetoelectric (ME) effect is associ- 
ated with two types of exchange interaction. In collinear 
spin structures the ME coupling is explained by changes 
of the isotropic excha nge interaction leading to the ex- 
change strictionpEGD whereas in spiral spin structures 
the antisymmetric part of the anisotropic exchange inter- 
action is held responsible.^^ Despite the generally ac- 
cepted phenomenological description} 3 -' 5 -' the microscopic 
picture of the ME mechanism is much more complex and 
still lacks a unified explanation. 

In this paper we fo cus on the mechanism of the ME 
coupling in FeT^OsBrj 12 ! 13 ! where electric polarization in 
the long-range ordered elliptical IC amplitude-modulated 
(AMOD) magnetic state was suggested to originate from 
phase shifts between the exchange-coupled AMOD mag- 
netic waves and thus argued to differ from conventional 
ME mechanisms. The system adopts a layered struc- 
ture of [Fe4Oi6] 20_ tetramer clusters connected via Te 4+ 
ionsP^l The magnetic lattice is composed of alternat- 
ing antiferromagnetic Fe 3+ (S = 5/2) spin chains cou- 
pled by frustrated interactions,^ which amount to ~ 1/3 
of the dominant intrachain interaction J2~19K. Two 
subsequent magnetic transitions were identified.^ The 
first, from paramagnetic to the high-temperature IC 
magnetic state (HT-IC) with a constant wave vector 
Qici = [\ 0.466 0), occurs at T/vi = llK and is rapidly 
followed at T^2 = 10.5 K by the second one into the low- 
temperature IC multiferroic state (LT-IC). The ellipti- 



cal IC AMOD order in the LT-IC phase is characterized 
by long axis of the ellipsis along the (1,-1, 0.2) direc- 
tion in the a* be orthonormal system (used throughout 
the paper) and the magnetic wave vector that pro- 
gressively changes from qici to qic2 = (5 0.463 0), where 
it settles below T ~ 6 KP^The accompanying electric po- 
larization, ascribed to the exchange striction of the inter- 
chain interactions, points along the c-axisP^ In contrast 
to the LT-IC state, the magnetic ordering in the HT-IC 
phase is still unknown and thus hampers the understand- 
ing why the electric polarization is absent in this phase 
and why it develops in the LT-IC phase™ 

Using combined nuclear magnetic and quadrupolar res- 
onance (NMR and NQR, respectively) , spherical neutron 
polarimetry (SNP), and neutron diffraction techniques, 
we solved the magnetic structure in the HT-IC phase. 
Here the magnetic moments are, like in the LT-IC phase, 
sinusoidally modulated, but are now almost completely 
collinear with the b axis. Furthermore, our NQR results 
indicate changes of the electric-field gradient (EFG) at 
the Br sites below T^2 (in the multiferroic LT-IC phase) , 
corroborating the minute displacements of the Tc 4+ ions 
that manifest as a bulk electric polarization.^ Compar- 
ison with the LT-IC magnetic structure implies that in 
the HT-IC state the phase shifts between certain mag- 
netic AMOD waves are suppressed in acco rdance with 
the proposed ME coupling mechanism! 12 ! 15 ! These phase 
shifts are thus most likely responsible for the lack of the 
electric polarization above T/V2- 



II. EXPERIMENTAL DETAILS 

NMR and NQR measurements were performed on 
high-quality single crystals^ with an average size of 
15x8x2 mm 3 on a home-build spectrometer in the tem- 
perature range between 4 and 300 K in zero magnetic 
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TABLE I. Upper panel: DFT calculated components of the 
electric-field gradient (EFG) tensor V%j given for the Bri 
and Br2 sites in the a* be coordinate system in units of 
10 21 V/m 2 . Lower panel: the corresponding quadrupole split- 
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field and at 4.7 and 9.4 T. 

Spherical neutron polarimetry (SNP) was performed 
at 10.7 K on the same crystals using a MuPAD device 
on the triple axis spectrometer TASP (A = 3.2 A) at the 
Swiss Neutron Spallation Source (SINQ), Paul Scherrer 
Institute (PSI), Switzerland. Intensities of the magnetic 
reflections at temperatures between 9.2 and 11 K were 
collected at the same location using the single-crystal 
diffractometer TriCS (A = 2.32 A). 



III. RESULTS 

A. NQR and NMR experiments 

Local-probe NQR and NMR experiments on 79 Br and 
81 Br nuclei with the spin J = 3/2 were chosen because, 
(i) in addition to a standard dipolar magnetic moment 
through which they detect the local magnetism, (ii) they 
also possess a quadrupole moment, making them sensi- 
tive to EFG and thus highly susceptible even to the tini- 
est structural deformations. From the experimental point 
of view, however, this makes it very difficult to find the 
resonant frequency, since EFG can vary within several 
orders of magnitude, depending on the details of the lo- 
cal Br environment. In addition, in FeT^OsBr there are 
two crystallographically inequivalent Br sites; Bri that 
is coupled to a single magnetic Fe 3+ (5 = 5/2) ion, and 
Br2, interacting with three magnetic Fe 3+ ions [inset to 
Fig-Bb)]. 



1. Characterization of local Br environment 

We first performed density functional theory (DFT) 
calculations of the EFG tensors at the Br! and Br 2 
sites (for details see Appendix[B|) , in order to facilitate 
the search of the 79 Br and 81 Br NQR and NMR sig- 
nals. The obtained EFG tensors (Table]!]) imply that the 
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is significantly larger than at the Br 2 site. Here, 79,81 



quadrupole splitting 



180 



150 



120 



90 
12 



60 



120 



180 240 



300 




(b) 






o 5 Fe ' 


B = 4.7T 








- v = 58.9 MHz As/ 










OA 

o Br , Fe - 


0>0 F e 2 . 









10 



12 14 
T(K) 



16 



20 



FIG. 1. (Color online) (a) Angular dependences of the sl Br 
and the 79 Br NMR central lines (— 1/2 o 1/2) at 80 K and 9 T 
(symbols) for the Bri site and simulations based on DFT cal- 
culated EFG tensor (lines). Inset: The corresponding NMR 
spectrum of both isotopes at </> aca = 120°. (b) Temperature 
dependence of the NMR relaxation rate 1/Ti at the peak of 
the central 81 Br 2 line, i.e., at 58.9 MHz, in the field of 4.7 T 
along a*. Solid lines are guides for the eyes. Inset: local 
coordinations of the Bri and Br 2 sites. 



denotes the quadrupole moment of the 79 > 81 Br isotopes, 
h is the Planck's constant, e is the electron charge, while 
Vij (i, j — x, y, z) are the components of the EFG tensor. 
The EFG tensor for Br! is almost axially symmetric with 
the asymmetry parameter 7] = {V XX — V yy )/V zz = 0.09, 
whereas for the Br 2 site 77 = 0.58. The 81 Br NQR sig- 
nals at 260 K (deep in the paramagnetic phase) were ex- 
perimentally found at 165.0 and 38.4 MHz for Br! an d 
Br 2 , respectively (Fig. [2]), which is within ~10% of the 
valued^ vq(1+t] 2 /3)s predicted by DFT calculations. 
The remarkable accuracy of the DFT results was further 
tested by measuring angular dependences of the 79 ' 81 Br 
NMR signals for Bri around all three crystallographic 
axes (a*, b, and c) in the field of 9 T at 80 K [Fig.Qa)] . If 
Bri EFG values are increased by 8.2 % with respect to the 
DFT calculations, the angular dependence, calculated by 
exact diagonalization of the nuclear spin Hamiltonian (for 
details see Appendix|A~|), nicely matches the experiment 
[lines in Fig.JTJa)]. We note that small discrepancy be- 
tween calculations and experimental data probably orig- 
inates from additional hyperfine fields and/or tiny mis- 
alignment of the crystal. 

Having determined the quadrupolar interactions, it is 
now our task to clarify if 81 Br NMR probes the mag- 
netism as well. We measured spin lattice relaxation 1/Ti, 
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FIG. 2. Temperature dependence of the 81 Br NQR signal 
for the Bri and the Br2 sites. Insets: Corresponding spectra 
measured at 80 K. 



FIG. 3. (Color online) Normalized 81 Br NQR spectra for (a) 
the Bri and (b) the Br2 sites. Different colors (shades) corre- 
spond to the three different magnetic phases (paramagnetic, 
HT-IC and LT-IC). 



which is a highly sensitive parameter for the critical spin 
fluctuations in the vicinity of the magnetic transitions. 
The temperature dependence of 1/Xi for the 81 Br2 cen- 
tral NMR transition [Fig.[]|b)], at 58.9 MHz, clearly re- 
flects two distinct lambda-type anomalies, signifying the 
two magnetic transitions. Compared to the zero field re- 
sults the splitting between the two transitions increases 
by ^0.5 K, due to the magnetic field of 4.7 T applied 
along the a* axisP^ 



2. Low temperature NQR spectra 

On cooling from 260 K both 81 Br NQR lines (Bri 
and Br2) shift linearly to higher frequencies down to 
~30K, where they are no longer temperature depen- 
dent (Fig.[2|. This shift is a result of a slight (-2%) 
increase of the EFG values due to the crystal lattice con- 
traction. As expected, below the first magnetic tran- 
sition, at T/vi = llK, the intensity of the sharp para- 
magnetic resonance gradually transfers to a broad U- 
shaped signal (Fig. [3]), typical for a sinusoidal distri- 
bution of local magnetic fields in the IC structures.^ 
Clearly, the two signals coexist in a narrow tempera- 
ture region around Tjvi, i.e., approximately between 11.0 
and 10.6 K, revealing the first-order nature of this tran- 
sition. To clarify the origin of the IC modulation in 
the HT-IC phase, i.e., whether is it solely magnetic or 
also structural, we in parallel measured NQR signals for 
the Br isotope, which has smaller gyromagnetic ratio 
7 and larger quarupolar moment Q than the 81 Br iso- 
tope ( 79 7 = 10.6663 MHz/T < 81 7 = 11.4978 MHz/T and 
79 Q = 31.3xl0- 30 m 2 > 81 Q = 26.2xl0- 30 m 2 ). Over- 
plotting the 79 Br and the 81 Br signals in the HT-IC phase 
we find that for both sites their widths scale with 7's 
[Figs.[4ja),(b)] and not with Q's. This proves that the 
observed IC modulation is solely magnetic. In addition, 



the center of gravity of the lines does not shift at T/vi 
(Fig. [3]), indicating that the EFG tensors are unaltered, 
i.e., suggesting that magnetic transition at Tjyi does not 
induce any significant crystal structure distortions. 

On further cooling, both signals dramatically change 
again at T/v 2 = 10.5 K, i.e., at the transition from the HT- 
IC to the LT-IC phase. Clearly, to achieve proper scaling 
of the 79 Br and 81 Br LT-IC signals with 7's [Fig.gd)], 
vq at the Br 2 site has to be reduced by ~2% [inset to 
Fig.Qd)] . This implies that EFG at the Br 2 site is sensi- 
tive to minute lattice distortions, accompanying the elec- 
tric polarization in the LT-IC phase. 12 ' Modification of 
vq for the Br! site is to small to be assessed from our 
measurements [Fig.Qc)]. In addition, as opposed to the 
HT-IC phase where both Bri and Br 2 sites have simple 
U-shaped NQR spectra, in the LT-IC phase their spectra 
suddenly become completely different (Fig. [3]). Further 
splitting of the Br 2 spectra below T^2 implies the loss of 
certain symmetries in the multiferroic state that might 
be still present in the HT-IC phase. 

B. Neutron diffraction experiments 

1. Magnetic structure in the HT-IC phase 

To verify whether the different NQR spectra originate 
from different magnetic-order symmetries in the HT- 
IC and the LT-IC phases, we next decided for neutron 
diffraction experiments, aiming to determine the mag- 
netic order in the HT phase. The combination of spher- 
ical neutron polarimetry (SNP) and conventional single 
crystal neutron diffraction has proven very useful in the 
past for determination of complex magnetic structures, 
e.g., IC arrangements or systems with superimposed nu- 
clear and magnetic contributions.^ Compared to conven- 
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FIG. 4. (Color online) The 81 Br and 79 Br NQR spectra for 
(a) the Bri and (b) the Br2 sites in the HT-IC phase and for 
(c) the Bri and (d) the Br2 sites in the LT-IC phase. In the 
inset a scaling with reduced v'q =0.98i/q is shown for Br2. 



TABLE II. Components of vectors Sg m = (Sq x , So y , Sq x ) for 
s = Re, Im, defining the elliptical envelops for two independent 
magnetic atoms (Fei and Fe2) for the best magnetic structure 
model at 10.7 K, and eight magnetic phases i/> mn in units of 
2-7T, i.e., one for each of the magnetic Fe mn atoms in the unit 
cell (m — 1,2, n= 1-4). The sites Fei2-Fei4 are obtained from 
Fen [0.1184(6), -0.001(1), 0.9734(7)] and Fe 22 -Fe 2 4 fromFe 2 i 
[0.9377(6), 0.2953(1), 0.8562(6)] by symmetry elements i, 2 ly 
and 2i y i, respectively. The orientation of the moments is 
given in the a*bc coordinate system, while |Srj| ~ 1.2 \xb- 
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tional single-crystal neutron diffraction, the SNP method 
has enhanced sensitivity to the direction of the magnetic 
moments and thus allows to differentiate between com- 
plex magnetic structures, e.g., between AMOD and he- 
lical spin arrangements. The two experiments were con- 
ducted at 10.7 K, with SNP performed for three different 
crystal orientations. In addition to the hkO orientation, 
where the scattering plane was defined by the (10 0) and 
(0 1 0) reciprocal vectors, the crystal was rotated to the 
scattering plane defined by the (0 1 0) and (1 2) vectors, 
and finally to the scattering plane with (0 1) and either 
the (0.5 0.466 0) or (0.5 0.534 0) vectors. Altogether, we 
accumulated 24 polarization matrices and 62 integrated 
intensities. 

Starting with the representation analysis, we find that 
magnetic wave vector qici = (h 0.466 0) breaks the inver- 
sion symmetry already in the HT-IC phase. This leaves 
two possible one-dimensional irreducible representations 
of the little (magnetic) group, which couple magnetic mo- 
ments at the Fe sites related by a 2i y twofold screw axisP^ 
Since the presence of the 2± y symmetry would explain the 
lack of the electric polarization (in the ac plane) as well 
as the high symmetry of the NQR spectra in the HT- 
IC phase, we start the refinement of the corresponding 
magnetic order assuming a single irreducible representa- 
tion. In the most general case, the magnetic moment at 
a particular Fe site is defined as 



S m nOi) = S^ n cos(q-r l --0 mn ) + So" n „sin(q-r i -V'mn). 

(1) 

Here, the vector r, defines the origin of the i-th unit 
cell, to = 1,2 identifies the crystallographically inequiv- 
alent Fe-sites, and 77=1-4 denotes the four Fe positions 
within the crystallographic unit cell (for details see the 
caption of Table H| . The complex vector Sq mn is deter- 
mined by its real and imaginary components, and 



S]^ mn , which define the amplitude and the orientation of 
the magnetic moments, i.e., the envelope of the magnetic 
cycloid/spiral, while ip mn denotes its phase shift. The 
magnetic wave vector q is in units of (2ir/a, 27r/6, 2ir/c) 
and ip mn in 27r. We stress that within a single ir- 
reducible representation S 0m („+2) — (=tl, Tl> ±1) ■ Somn 
and ^m(n+2) = i>mn + q y lci /2 1 where 77 = 1,2. To avoid 
overparameterization of the problem we assume the same 
complex vector Srjmi = Srjm2 = So m - Moreover, to assure 
the best assessment of the experimental uncertainty, the 
estimated standard deviations of the polarization matri- 
ces and the overall refinement were treated in the same 
way as in our study of the LT-IC phaseP^ Surprisingly, 
neither of the two irreducible representations can de- 
scribe the HT-IC data satisfactory, as both refinements 
diverge, implying that all symmetry operations are bro- 
ken already in the HT-IC phase. 

In the next step we, therefore, resort to the elliptical IC 
structure model used to describe the LT-IC phase. Here, 
the symmetry relations between S 0m ( n +2) an d Sg ra „ 
and between ip m (n+2) and ipmn f° r 77 = 1,2 do not ex- 
ist anymore, while we extend the relation So m n = So m 
to include 77 = 3,4. In addition, as in the LT-IC casej^ 
we allow different domain populations for each experi- 
ment. Indeed, this refinement leads to a stable solution, 
which is almost completely sinusoidally modulated with 
Sq 1 ^ 3> Sq£i and the dominant components (Sg™J of the 
magnetic moments aligned very close to the b axis (Fig. [5] 
Table[H]) . Goodness of the refinement reflects in the to- 
tal cost Ctot =J2j Xj / ( N j obs-N par )= 46.6, which is close 
to C to t = 33.7 obtained in the refinement of the LT-IC 
structure^! and better than C to t = 57.3 for a simplified 
sinusoidal (collinear) AMOD model of the HT-IC phase. 
It is significantly better than Ctot = 169.3 for an alter- 
native circular cycloidal model. In the above expression 
for Ctot j is the number of the datasets and A^par is the 




FIG. 5. (Color online) Evolution of Fen and Fe2i magnetic moments in the a*b projection along the b axis in the HT-IC [light 
(red) arrows] and LT-IC [dark (blue) arrows] phases. The rectangle represents the unit cell. Note that the elliptical envelop 
for the HT-IC phase is almost completely flattened, with moments pointing approximately along the b axis. 



number of the fitting parameters. For each dataset with 

iV iobs observations Xj =Ei^i° b3 (^ob s -^caic) 2 / cr i iob3 : 
where X denotes polarization matrix elements P or inte- 
grated intensities /, and <Jx iabB is the estimated standard 
deviation of the observation. We point out that the re- 
finement of the LT-IC phase included larger Nj ^ s , while 
iVp ar was the same, which led to a somewhat lower Ctot- 
In short, the results of our refinement surprisingly show 
that magnetic ordering removes all symmetry restrictions 
for the electric polarization already in the HT-IC phase. 



2. Temperature 



In order to understand why there is no electric po- 
larization in the HT-IC phase, we next measured tem- 
perature dependence of 20 magnetic reflections between 
11.0 K and 9.2 K. This allows us to follow the evolution 
of the magnetic structure during the transition from the 
HT-IC to the LT-IC phase. Due to the limited amount 
of data, rather small ordered magnetic moments in the 
investigated temperature interval and since the simpler 
collinear AMOD model already captures the most essen- 
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FIG. 6. (Color online) Temperature dependence of the M x 
and My components of the magnetic moments deduced from 
the intensities of the neutron reflections on heating through 
the transition from the LT-IC to the HT-IC phase. Solid lines 
are guides for the eyes. 



tial properties of the more complicated elliptical model 
our refinement considers the collinear sinusoidal AMOD 
magnetic structure model. In particular, we assume that 
magnetic moments lie in the a* b plane and are strictly 
AMOD with the same S = S^f 1 for all Fe sites. This 
way we focus on magnetic phase shifts ip mn , which were 
pointed out as potential source of the ME coupling in 
earlier studies.^ The results clearly show that with in- 
creasing temperature the a* component (M x ) is reduced 
and completely disappears at Tjv2, whereas the b compo- 
nent (My) starts to decrease only in the vicinity of T/vi 
(FigM. This indicates that the Fe 3+ magnetic moments 
turn from the (1-10) direction in the LT-IC phase to- 
wards the (0 1 0) direction in the HT-IC phase. On the 
other hand, the changes of V'mn are minute and appear to 
be insensitive to the LT-IC to HT-IC transition. In fact, 
reducing the number of independent ip mn only mildly af- 
fects the quality of the refinement, i.e., R factor increases 
from ^6 to ~9, and it does not affect the derived rota- 
tion of the magnetic moments. This suggests that our 
temperature-dependent-neutron diffraction data are in- 
sufficient to reliably extract the temperature evolution of 
the magnetic phases. Similarly, we were unable to refine 
the temperature dependence of the tiny magnetic com- 
ponent along c (M z ). 



IV. DISCUSSION 

A. Long-range magnetic ordering 

At first sight, the results of the two complementary 
experimental techniques are contradictory. The neutron 
diffraction suggests that magnetic ordering breaks all 
crystal-symmetry relations already in the HT-IC phase, 
whereas NQR implies that the HT-IC phase is more sym- 
metric than the LT-IC phase. In order to clarify this issue 
and to extract as much information about the long-range 
ordering in FeTe205Br as possible, hyperfine coupling 
tensors have to be determined. Since both Br sites (Bri 
and Br 2 ) lie at general positions and since the hyperfine 
coupling interaction is symmetric in the first order, we 
need to find for each hyperfine coupling tensor all six 
components. For simplicity we start with Bri, which is 
coupled to a single Fe ion and thus its hyperfine inter- 
action depends only on one hyperfine coupling tensor. 
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TABLE III. Derived hyperfine coupling tensors for the Bri 
and Br2 sites in the a* be coordinate system in units of mT//xs. 
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FIG. 7. (Color online) The 81 Br NQR spectra for (a) Bn 
and (b) for both isotopes of Br 2 in the HT-IC and the LT-IC 
phases. Thick black lines are simulations considering mag- 
netic structures determined by neutron diffraction and de- 
rived hyperfine coupling tensors (see text for details). For 
comparison we show simulations of the HT-IC spectra based 
on the LT-IC structure with appropriately reduced size of the 
magnetic moments [thick light (magenta) line]. 



Taking into account the LT-IC magnetic structure^ and 
the obtained EFG tensor (Figs.[l][2]), we determine the 
Br! hyperfine coupling tensor by fitting the NQR Brx 
spectrum measured at 4 K (for details of the NQR spec- 
trum calculations see Appendix A}. We stress that the 
derivation of the hyperfine coupling tensor from the an- 
gular dependences of the paramagnetic 79 > 81 Br NMR sig- 
nal [Fig.jlja)] has been avoided due to possible crystal 
misalignments (<5°), which can for so large EFG's re- 
sult in resonance shifts of several MHz. The obtained 
Bri hyperfine coupling tensor (Table [lH"| ) yields a good 
agreement between the experimental and calculated spec- 
tra [bottom panel in Fig.[7ja)]. Moreover, considering 
the HT-IC magnetic structure (Table|ll|) and the derived 
hyperfine coupling tensor (Table|H~ij), the HT-IC 81 Bri 
NQR spectrum is reproduced with a high accuracy [top 
panel in Fig.^a)] with no adjustable parameters. This 
unambiguously validates the orientation and amplitude 
modulation of the Fe 3+ moments in the HT-IC magnetic 
structure, as determined by neutron diffraction. Quite 
importantly, if the small c (M z ) component of the mag- 
netic moments is neglected, the HT-IC spectrum cannot 
be reproduced satisfactory, whereas disregarding the M x 
component (along the a* ) or ellipticity has almost no ef- 
fect on the simulated spectrum. 

Encouraged by a very good agreement between the 
neutron scattering results and NQR data for the Br! 
site, we focus now on the more complicated Br 2 NQR 
spectrum. Applying the same procedure, i.e., fitting of 
the LT-IC NQR spectrum by adjusting the hyperfine 
coupling tensor (the EFG is reduced by 14 % compared 
to DFT calculations) and considering the coupling with 
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three different Fe 3+ sites, each with its own hyperfine- 
coupling tensor, we manage to reproduce the main fea- 
tures of the Br 2 NQR spectrum as well [bottom panel in 
Fig.^b)]. To ensure a proper scaling of the EFG tensor, 
spectra for both isotopes were fitted simultaneously. Our 
fitting results get even greater value when one considers 
that during such a broad frequency sweep the measured 
NQR intensity may significantly vary due to frequency- 
dependent sensitivity of the spectrometer and the use of 
several experimental setups. In addition, we note that 
the discrepancy between the experimental and the calcu- 
lated spectrum can be also due to tiny modulation of the 
vq, which can result from the weak IC structural mod- 
ulation found in SNP study of the LT-IC phaseP This 
further complication is beyond the scope of our simula- 
tions, but may explain why the 79 Br2 and 81 Br2 NQR 
spectra scaled by 7's do not match perfectly even when 
reduced vq in the LT-IC phase is considered [inset in 
Fig.[4jb)]. Nevertheless, the derived hyperfine coupling 



tensors (Table III ) can reproduce the width and the main 
spectral singularities also for the HT-IC phase [black lines 
in the top panel in Fig.[7^b)] when the EFG is increased 
by 2 % compared to the LT-IC phase, in accordance with 
our previous observations [Figs.|4][b),(d)]. 

For comparison, we plot also calculated spectra corre- 
sponding to the LT-IC magnetic structure with appro- 
priately scaled magnetic moments to match the HT-IC 
values. For these the discrepancy from the experimental 
data is much more pronounced for both Br sites [light 
lines in Fig. [7]. The overall agreement between the neu- 
tron diffraction and the NQR results thus offers a confir- 
mation of the proposed magnetic structures in the LT-IC 
and the HT-IC phases. It also shows that the simplicity 
of the Br 2 NQR spectra in the HT-IC phase originates 
from a smaller size of the ordered magnetic moments and 
a changed orientation of spins, and not from the higher 
symmetry of the magnetic structure. Finally, we note 
that the similarity between ip mn and VVi(n+i) f° r ^=1,3, 
i.e., between the sites, related by the inversion symmetry, 
implies that even though the inversion symmetry is bro- 
ken already by qici, the system effectively reduces this 
effect by matching the relevant phases. 
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B. Magnetoelectric coupling 

The most important experimental finding of our study 
is that all crystal-symmetry relations are broken already 
by the HT-IC magnetic order, which should thus, in prin- 
ciple, also allow for the establishment of the electric po- 
larization. This is in line with the first-order nature of 
the transition, suggested by the coexistence of the HT- 
IC and the paramagnetic phase in a narrow tempera- 
ture range around Tjvi (Fig. [3]). Surprisingly, the elec- 
tric polarization does not develop until the second mag- 
netic transition into the LT-IC phase. In addition, the 
lack of the symmetry relations removes all limitations re- 
garding its orientation, leaving no clue why the actual 
polarization points along the c axis. The observed re- 
sponse, therefore, deviates from other multiferroics where 
the multiferroic phase evolves in two subsequent [or one 
as in RbFe(MoG'4)2p' continuous magnetic transitions 
and thus allows to exploit the phenomenological descrip- 
tion of the ME coupling to its full extent In particu- 
lar, in contrast to our case, continuous transitions pre- 
serve certain relations between the magnetic ordering in 
the multiferroic phase and the symmetries of the crystal- 
lographic space group, which enables predictions of the 
direction of the emergent electric polarization!^ 

In FeTe205Br, the electric polarization could, in prin- 
ciple, be associated with the reorientation of the magnetic 
moments from the (1-10) direction in the LT-IC phase 
towards the (0 1 0) direction in the HT-IC phase. How- 
ever, such scenario would be most probably associated 
with the simultaneous reorientation of the electric polar- 
ization, which contradicts our previous results (Ref. ll2p 
showing that the orientation of the electric polarization 
is temperature independent in the entire LT-IC phase. A 
second possibility is that electric polarization is associ- 
ated with the ellipticity in the LT-IC magnetic structure, 
which is again not very likely, as our refinement shows fi- 
nite ellipticity also in the HT-IC phase. 

To understand the absence of the electric polarization 
in the HT-IC phase, we thus compare the LT-IC and 
HT-IC magnetic structures in respect to the so-called 
magnetic-phase-shift ME coupling mechanism, which 
predicts electric polarization PocMi ■ Mj sin(Aipk)^ 
Here Mij are the magnetic order parameters correspond- 
ing to a pair of the exchange-coupled AMOD magnetic 
waves, while Aipk denotes the phase shift between the 
two. Generally, each of the six possible exchange inter- 
action J k (k = 1-6^ can be involved in the ME coupling 
mechanism. We find, however, that on heating from the 
LT-IC to the HT-IC state Aipk changes towards ir or 
to 0, i.e., leading to P—^0, only for k — 4. Similar, but 
significantly weaker, trend is noticed for J5, which has 
been also highlighted before as the most likely candidate 
to drive the exchange strictionP^l On the other hand, all 
changes corresponding to other fc's are rather irregular. 
These observations are further supported by the fact that 
the NQR spectrum for Br 2 , which is coupled to the J4- 
bridging Te4 ion, shows a pronounced change of vq and 



thus reveals tiny structural transformations, most likely 
related to the onset of the electric polarization. 

The above argumentation, therefore, suggests that the 
most coherent explanation of the magnetic ordering and 
its relation to the ME coupling is provided by the "mag- 
netic phase shift" mechanism involving J4 and possibly 
J5 exchange pathways. 



V. CONCLUSIONS 

We have investigated the magnetic ordering in the HT- 
IC phase of the FeTe2 05Br system by combining com- 
plementary neutron diffraction, nuclear quadrupolar and 
magnetic resonance techniques. We find that due to 
the first-order transition from the paramagnetic phase 
all crystal symmetries are broken already in the HT- 
IC phase, which makes it, from the symmetry point of 
view, equivalent to the multiferroic LT-IC phase. How- 
ever, the ellipticity in the HT-IC phase is significantly 
reduced and the magnetic moments with sinusoidally 
modulated amplitudes align almost exactly along the b 
axis. Furthermore, the phase shifts between the magnetic 
AMOD waves, corresponding to the J4 exchange interac- 
tion, converge towards Aip^ — > or 7r, which according 
to the magnetic-phase-shift ME coupling mechanism^ 
yields P cx Mi ■ M ; sin(A?/>fc) — > 0, and could thus ex- 
plain why the electric polarization vanishes in the HT-IC 
phased In addition, changes of the EFG at the Br 2 site 
in the LT-IC phase imply minute displacements of the 
Te 4+ ions, which must be associated with the emergent 
electric polarization. 
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Appendix A: Calculations of NMR and NQR spectra 

Usually the NMR and NQR absorption lines are cal- 
culated for the two limiting cases. In case of NMR ex- 
periment, Zeeman term in the nuclear spin Hamiltonian 
is typically taken as the dominant one, whereas electric 
quadrupole effects are considered as a perturbation. On 
the other hand, for NQR, so-called Zeeman perturbed 
electric quadrupole Hamiltonian is assumed. In the case 
of FeTe2 0sBr, however, the NQR frequencies are com- 
parable to the NMR ones and thus exclude the possi- 
bility to use simple perturbative approaches. Therefore, 
our approach to calculate the resonance frequencies and 
their intensities is based on the exact diagonalizatiorP2l 



of the complete Hamiltonian for magnetic resonance of 
quadrupolar nuclei 2 ^ 

H = Hz + H Q +H hyp + H di p. (Al) 

Here Hz denotes the Zeeman term, Hq the quadrupole 
interaction, and "Hh yp and "Hdip the influence of the hy- 
pcrfinc and dipolar fields, respectively. The individual 
terms have the following form, 

Hz = -lhB -l, (A2) 

H Q = 4J(2 j _ 1) [ y o(3^ ~ I 2 ) + V +1 (I-I Z + IJ-) 

+V. 1 (I+I Z + IJ+) + V +2 (r) 2 + V. 2 {I + ) 2 }, 

(A3) 

H hyp = -^H{S)-A-I, (A4) 

W di p = -7fiB dip .I. (A5) 

Here 7 denotes the nuclear gyromagnetic ratio, 
I=(I x ,I y I z ) is the nuclear spin, I ± =I x ±iI y , A is 
the hyperfine coupling tensor, Vb = V zz , V±% = V zx ±iV zy , 
V±2 = h(V xx — Vyy) ±iV xy with Vij being the components 
of the EFG tensor, (S) is the time averaged electron 
magnetic moment, and B^ip and Bo are the dipolar and 
the applied external magnetic fields, respectively. In the 
long-range ordered magnetic states (S) and Bdi p are ex- 
actly determined by the LT- and HT-IC magnetic struc- 
tures given in Ref.[T3]and Table|Tlj respectively, whereas 
in the paramagnetic phase magnetic moments are as- 
sumed to lie along Bo with amplitudes scaled by the 
magnetic susceptibility.^ Exact diagonalization of the 
above Hamiltonian allows calculation of the nuclear spin 
eigenstates and thus enables to extract corresponding 
NMR/NQR transition frequencies for any kind of lo- 
cal (external or internal) magnetic field and for an ar- 
bitrary EFG tensor. In addition, considering the orien- 
tation of the excitation/pick-up coil one can derive also 
the probabilities of the individual magnetic transitions 
and can thus estimate the intensities of the correspond- 
ing NMR/NQR absorption lines. By rotation of the crys- 
tal system with respect to the laboratory system, the 
NMR/NQR spectrum, which for Br nucleus with 1 = 3/2 
in most general case consists of six absorption lines, can 
be calculated for any orientation of the crystal. 

The 79 - 81 Br NQR spectra in the IC long-range ordered 
magnetic phases were computed for each of the two mag- 
netic domains (NMR/NQR experiments differentiate the 
two domains related by 2± y symmetry) by summing the 
four (n = 1-4) spectral contributions for one crystal unit 
cell and than by summing the contributions of 100 con- 
secutive cells along the IC direction (b axis), as shown for 
a single domain in Fig. [8] for 81 Br. The complete spec- 
tra (Fig. [7]) were finally obtained as a sum of two contri- 
butions corresponding to the two equally populated do- 
mains. We note that for each crystallographicaly unique 
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FIG. 8. (a) The simulated spectra for individual 81 Br2 n sites 
for a single domain, (b) Corresponding distributions of the 
amplitude of the local magnetic fields Bi oc = — 7/i((S) ■ A + 
Bdip) for the four contributing transitions between different 
nuclear spin states. We note that the intensity of the calcu- 
lated spectra reflect the density distribution of the local fields, 
i.e., the peaks in the spectra coincide with the most common 
local fields. 

Br site (Br! or Br 2 ) electronic dipolar fields were calcu- 
lated individually at each Br nucleus considered in the 
above summation (2x4x 100) by assuming a sphere large 
enough (~40A) to ensure convergence. 



Appendix B: DFT calculations of EFG tensors 

The components of the EFG tensor were calculated ab 
initio within the framework of the density- functional the- 
ory by applying the Wien97 codep^ which adopts the full- 
potential linearized-augmented-plane-waves (FLAPW) 
method: 26 The experimental data for the lattice param- 
eters and the atomic positions served to describe the in- 
put crystal structure at room temperature,^ whereas the 
muffin-tin radia were 2.1a.u. for the Fe atoms, 1.98a.u. for 
the Te atoms, 1.5a.u. for the O atoms, and 2.68a.u. for the 
Br atoms. The exchange-correlation effects were treated 
within the local-density approximation (LDA)f2^The in- 
tegration over the Brillouin zone (BZ) was discretized 
by summing up 333 k-vectors in terms of the Gaussian 
methocp3! with the smearing parameter of 0.02 Ry. The 
plane- wave-expansion cut-off energy was set to 16 Ry and 
the magnitude of the largest wave vector in the Fourier 
expansion of the charge density was 10a.u. _1 . The EFG- 
tensor components are defined as the second derivates 
of the Coulomb potential at the particular nucleus. The 
Coulomb potential is obtained from the total charge den- 
sity by Solving the Poisson's equation. The calculation 
of the EFG is therefore straightforward once the non- 
spherical components of the charge density p(r) are avail- 
able as it is the case in the FLAPW method. 
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